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Abstract 

We consider a 6D space-time which is periodic in one of the extra dimensions and com- 
pact in the other. The periodic direction is defined by two 4-brane boundaries. Both static 
and non-static exact solutions, in which the internal spacetime has constant radius of cur- 
vature, are derived. In the case of static solutions, the brane tensions must be tuned as in 
the 5D Randall-Sundrum model, however, no additional fine-tuning is necessary between 
the brane tensions and the bulk cosmological constant. By further relaxing the sole fine- 
tuning of the model, we derive non-static solutions, describing de Sitter or Anti de Sitter 
4D spacetimes, that allow for the fixing of the inter-brane distance and the accommoda- 
tion of pairs of positive-negative and positive-positive tension branes. Finally, we consider 
the stability of the radion field in these configurations by employing small, time-dependent 
perturbations around the background solutions. In analogy with results drawn in 5 dimen- 
sions, the solutions describing a de Sitter 4D spacetime turn out to be unstable while those 
describing an Anti de Sitter geometry are shown to be stable. 
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1 Introduction 



It would be an understatement to say that the possibility of resolving the hierarchy problem 
in models with a warped extra dimension |1J has received considerable attention over the 
last two years. Indeed, brane world models have dominated the literature in high energy 
theory The reason lies in its simplicity. By postulating the existence of two 3-branes with 
non-zero tensions, separated along the extra dimension by a distance L, in the background 
of a non-zero (negative) cosmological constant, one finds a simple solution for the scale 
factor along the extra dimension, a(y), which is exponential. Thus, length scales (and 
hence mass scales) on one brane are exponentially enhanced (or suppressed) relative to the 
other. A mass hierarchy naturally arises between the two branes which can be labelled the 
Planck and weak branes respectively. 

Of course, there is a price to pay for this simplicity. First, as it is well known, the 
tensions of the two 3-branes, must be fine-tuned so that Ai = — A2. Second, these tensions 
must be tuned to the bulk cosmological constant, A B , in order to produce a static solution. 
The origins of these fine-tunings come about when one considers static solutions to the 5D 
equations of motion. The scale factor in the extra direction takes the form a(y) = e~ ky and 
the equations of motion require k 2 = —k 2 Ab/Q, where k\ is the 5D Newton's constant. We 
are, therefore, led to an Anti de Sitter 5D spacetime with A# < 0. By putting branes in 
the theory, and requiring that the warp factor exhibits periodic behavior along the extra 
dimension, we obtain the so called jump conditions which give [a r (y))i/ai = — ftgAj/3, where 
[a'] represents the difference in a' on the two sides of the brane. For one brane say with 
positive tension placed at the origin (y = 0), we see that k = K5A1/6, and for the second 
brane, located at an undetermined distance L, one finds k = — K5A2/6. Thus, we arrive at 
the conditions Ai = — A2 = y^As/z^f. Finally, the distance L is chosen so as to resolve 
the hierarchy problem by noting that masses scale as a(y). For other recent attempts at 
solving the hierarchy problem with extra dimensions see 0, BJ. 

Several extra-dimensional attempts at resolving the hierarchy, or the cosmological con- 
stant, problem have considered six- or higher-dimensional models ||, ||. Space-times with 
more than 1 extra dimension can allow for solutions with most appealing features, particu- 
larly in spacetimes where the curvature of the internal space is non-zero. These solutions, 
exhibiting either spherical or cylindrical symmetry with respect to the extra coordinates, 
can accommodate an exponential dependence on one of the extra coordinates, thus, re- 
sembling the 5D RS mechanism for the resolution of the hierarchy problem. In addition, 
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it turns out that such space-times can play an important role in relaxing the degree of 
fine-tuning in the RS models ||. Finally, these models can provide a framework in the 
context of which the stabilization of the radion field naturally takes place: for example, in 
Ref. ||, it was shown that in space-times with a constant spatial curvature of the internal 
dimensions, one can find solutions with a global minimum in the effective theory for the 
radion field. 

In this paper, we look for solutions to the 6D equations of motion based on an internal 
space of constant curvature. We first present an exact static solution where the warp factor 
depends on both extra coordinates and, hence, does not exhibit any spherical or cylindrical 
symmetry. The dependence on one of the two extra coordinates is a purely exponential one 
thus resembling the profile of the warp factor in the case of the 5D RS model. In analogy 
with RSI, the space-time contains two 4-brane boundaries with equal and opposite brane 
tensions. This configuration requires the same fine-tuning that exists in RSI model due to 
the jump conditions imposed at the boundaries. However, the solution does not contain the 
additional fine-tuning between the brane tensions and the bulk cosmological constant which 
is replaced by the fixing of the size of the extra dimension along which the 4-branes extend. 
The inter-brane distance along the remaining extra dimension remains arbitrary and it may 
be fixed only through the introduction of an additional mechanism for the stabilization of 
the radion field 0, ||. 

We then proceed to derive non-static solutions in the context of the same model. In 
this case, the exponential behaviour along one of the two extra coordinates changes to cosh 
or sink-like allowing for the accommodation of pairs of branes with positive tensions or 
positive-negative, respectively. The jump conditions lead to the fixing of the locations of 
the two branes along the same dimension and the fine-tuning between the brane tensions 
disappears rendering this solution totally free of any fine-tuning. 

Both of the above solutions, static or non-static, have been derived under the assumption 
that the extra space-time remains static. We formulate an "extremization" constraint that 
may serve as a consistency check for any 6D solution with a constant or non-constant radion 
field. We finally perform a stability analysis around the solutions with a constant "extra" 
scale factor in order to check their stability under small time-dependent perturbations. We 
find, in agreement with similar results derived for 5D spacetimes || - ||14|| , that the system 
of two Minkowski 4D subspaces has a vanishing radion mass, a pair of two de Sitter ones 
has a negative mass squared while the system of two Anti de Sitter 4D subspaces has a 
positive mass squared. 
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In the next section, we present the model and derive the exact static solution in section 3. 
We show explicitly how the correlation between the brane tensions and the bulk cosmological 
constant is replaced by the fixing of the size of one of the two extra dimensions. In section 

4, we show how the relaxation of the fine-tuning between the brane tensions leads to de 
Sitter or Anti de Sitter expansions in the 4D space-time in analogy with 5D models. The 
"extremization" constraint and the stability analysis of our solutions are discussed in section 

5. Finally, we present our conclusions in section 6. 

2 The general framework 

Let us start by presenting the theoretical framework and geometrical set-up of our model. 
We first write down the action, that describes the gravitational theory of a 6-dimensional 
spacetime filled with a bulk cosmological constant, as[] 



where n\ = 8ti/Mq. The line-element of the 6-dimensional spacetime is assumed to be of 
the form 



In the above, x^ and (6, <ft) denote the coordinates along the usual four and two extra, 
all (initially) non-compact, dimensions, respectively. The function a(t, 8, if) represents the 
warp factor multiplying the 4-dimensional line-element, b(t) denotes the scale factor that 
determines the dynamics of the 2-dimensional extra spacetime while f(9) parameterizes its 
internal curvature. For reasons that will become clear shortly, we will often refer to the 
^-dimension as the "longitudinal" one and to the (^-dimension as the "transverse" one. 

In this paper, we will focus on the determination of static and non-static solutions under 

the assumption that the scale factor along the extra dimensions remains always constant. 

As noted in the Introduction, we are not assuming any specific mechanism for radion 

stabilization which would generate additional stress-energy terms 0, ||. We will instead 

assume that the internal curvature of the 2-dimensional extra spacetime will naturally lead 

the system to solutions with a constant radion field that correspond to a minimum of 

the radion effective potential as in Ref. ||. The outcome of this attempt is not however 

straightforward: in ||, it was assumed that a and b are functions of the coordinates x M 

throughout this paper, we follow Wald's conventions |L5]]: The metric signature is tjmn = ( — i +) — j +) 
and the Riemann tensor is defined as R" P11V — d^T" — d v T° + ... . 





ds 2 = a 2 (t, 6, ip) rj^ dx»dx v + b 2 (t) [d9 2 + f 2 (9) dip 2 } . 
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only. None of the scale factors depend on the internal coordinates 9 and <p and, thus, their 
results are not directly applicable to any warped brane model. 

Upon variation with respect to the 6D metric tensor, the above action and metric ansatz 
leads to the following Einstein's equations in the bulk 

IP Qba 3d 2 a 2 f" 3 /9ax2 3a d 2 a 
00 ~b 2+ ~ba~ + ~a T ~¥J~ Wf 2 v^V ~ Wf 2 d^ 2 

3a f da 3 /da\2 3a d 2 a 2 2 A (0 Q \ 
~Wd9~V^d9~) ~WW 2=Ha B ' [ } 

_^_26_26a a^_2a f_ _3_/9a\ 2 3a d 2 a 
u ~ ~T~ ~ba~ + ~a^ ~ ~a~ + ¥ J + Wp^d^) + ¥f 2 d^ 2 

3af da 3 /da\2 3a d 2 a 2 2a 
+ Wd9 + »W + VW 2= - K ' akB > (2 ' 4) 

Q bb 2bba _3b 2 a ^ 6 ( da \ 2 + 4 d 2 a 

55 a 2 a 3 a 3 a 2 f 2 ^dp' af 2 dip 2 



4/' 9a 6 /9a\ 2 „ „ 

G ft/ 2 & 2fe/ 2 ba 36 2 / 2 d | 6 sday | 6/ 2 /9ax2 

66 a 2 a 3 a 3 a 2 ^dp' a 2 ^d9* 

46 9a 3d da 3 d 2 a 
Go5 = ^bd9 + ^d9~adm =0 ' (2 - 7) 

46 9a 3d 9a 3 <9 2 a 
ao ay? a 2 ay? a atcty? 

G „ = i£* i *° 0. (2.9) 
a j dip a d9dip 

As we mentioned above, both extra dimensions are initially non-compact. Here, we 
choose to compactify along the ^-coordinate by introducing two 4-branes at the points 
ip — ipi and ip — ip2 (and impose periodic boundary conditions). The brane energy- 
momentum tensors can be written as 

= % 6 ~ yi) diag(-A 4 , -A,, -A,, -A,, — Aj, o) , (2.10) 
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where % — 1, 2 denotes the two branes. Note that the requirement of satisfying 6D energy- 
momentum conservation, T^.j^ = 0, together with the jump conditions (discussed shortly) 
forces one to introduce 1) an "inhomogeneity" reflected in the prefactor of eq. fl2.10|) and 2) 
an "anisotropy" such that T 5 5( ^ ^ j>*W f or — q, 3. On the other hand, as expected, the 
energy-momentum tensor in the bulk is dominated by the presence of a smoothly distributed 
bulk cosmological constant, i.e. 

T^ {B) = diag(-A B , -A B , -A B , -A B , -A b , -A fl ) , (2.11) 

whose sign remains arbitrary at this point. Note that, in the absence of a radion potential, 
we do not include an additional contribution to the (55) or (66)-component here. 

Under the assumption that the two 4-branes are infinitely thin, the discontinuity in the 
first derivative of the metric tensor along the ^-coordinate creates a 5-function contribution 



to its second derivative. The jump conditions that follow by matching the coefficients 



of the ^-functions on both sides of Einstein's equations, provide constraints on the first 
derivatives of the warp factor at the location of the branes. One may easily conclude, 
from the absence of any delta-function source on the right-hand side of the (66)-component 
of Einstein's equations, eq. fl2.6f) , that there is no discontinuity in the first derivative of 
the warp factor along the ^-coordinate due to the smooth distribution of energy along 
the longitudinal extra dimension. Then, from eqs. (|2.3j ) and ( |2.5|) , which contain second 
derivatives with respect to ip, we obtain respectively the constraints 



[d 9 a] 



ba 



«6 A 



<p=ip^ 3 b a 



K 6 ~ 



ip=(Pi 4 



— 7 Aj . (2.12) 



From the above conditions, it is obvious that the (ii) and (55)-components of the energy- 
momentum tensor of each brane must satisfy the relation: A, = 4Aj/3. Thus, the brane 
tension along the "extra" ^-coordinate is clearly distinguished from the remaining four, as 
pointed out below eq. ( |2.10| ). Similar "anisotropic" brane tensions can be found elsewhere 



in higher dimensional models. 



As noted above, the inhomogeneous prefactor in eq. (|2.10|) comes from the energy- 
momentum conservation constraint on Let us assume for the moment that T N is 
of the form (g(9) 5(ip — <Pi)/b) diag(— A i; — A i; — A i5 — A i: — Aj, 0^. We have allowed g to be 
arbitrary and we have retained the common normalization (1/6). The N = conservation 
equation is identically zero for all solutions with 6 = 0, which includes all of the solutions we 
consider. Note that this constraint on 6 stems from the anisotropy Aj ^ Aj. The N = 1, 2, 3 
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equations are trivial. The N = 5 equation reduces to 

-(</ + 4- + 4)Ai + 4-Ai = 0, (2.13) 
a j a 

with all derivatives taken with respect to 9. With the anisotropy relation derived from the 
jump conditions this becomes 

(g'+- + ^)A t = 0. (2.14) 
a j 

As we will show in the next section, our solution further requires ^ = y , from which it can 
be easily seen that g = l/f 2 . We also note that without this tension prefactor the jump 
conditions would contain an explicit factor of f(9) necessitating a constant / and hence a 
flat internal space. 



3 Static 6D solutions 



In this section, we focus on the derivation of an exact 6-dimensional static solution with 
an exponential warp factor and a constant scale factor, b(t) = bo. In this case, the non- 
vanishing components of Einstein's equations (after dropping all terms with time deriva- 
tives) take the form 

f" 3 ,da,2 3 d 2 a 3/' da 3 ( das2 3 &a _ 
f + cPP \d<p) + af 2 dp 2 + af d9 + a 2 W + a d9 2 ~ ^V^, (3.1j 

6 fda^2 4 d 2 a Af da 6 ( da^2 

6 /(9a\2 6 /9a\2 A d 2 a 2 7 2 . ,„ . 



a 2 / 2 ^d<p' a 2 ^06' a d9 2 

Af da A d 2 a 
af dtp a 89 dip 



0. (3.4) 



We assume that the dependence of the warp factor on the two extra coordinates can be 
written in a factorized form: a(9,<p) = Q{9) Eq. (|3.4j) , then, provides a very strong 

constraint on the functions 0(0) and f{9) leading to the relation 9(6*) = A e f(9), where A g 
is a constant. The difference of eqs. ( |3.2|) and ( [3.3|) , in turn, leads to the result 



% = ff"-f = u\ (3.5) 
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where u is again an arbitrary constant. The above equation allows us to write the gen- 
eral solution for $(0) as a linear combination of increasing and decreasing exponentials. 
However, eq. ( |3.3|) , which may be written as 

Q (^) 2 + 6 f' 2 + 4ff" = -4blf 2 AB, (3.6) 

restricts the form of the solution by demanding to be a constant as well. As a result, 
we may write the solution for <3> in the form 

^(p)=A v e h ^, (3.7) 

where an integration constant. Eq. Q3.7| ), together with eqs. ( |3.5| )-( |3~fi| ), brings the 
equation for the remaining unknown metric function, f(6), to the form 

/' 2 - A 2 / 2 + u 2 = , (3.8) 

whose general solution can be written as 

= _L [ e ±A(«-«o) + ^^-ft,)] , (3.9) 
4A 

where 

A 2 = -^A B . (3.10) 

Note that the parameter A can take real or imaginary values for a negative or positive, 
respectively, bulk cosmological constant. We will comment on these two different options 
shortly. 

We can furthermore easily check that the combination of eqs. (|3.5p -( ^6[) trivially sat- 
isfies the remaining equation ( |3.1| ). We may, therefore, write the full solution for the warp 
factor, that multiplies the line-element of the 4-dimensional space-time, as 

a(G, <p) = ^r ^ [e ±x ^ + 4u; 2 e^-*°)] , (3.11) 
4A 

where a stands for the product A^Aq. For simplicity, we may choose to place one 4- 
brane at the point tp = and the second one at ip = L. Then, by imposing the condition 
a(9 = 9q, if = 0) = 1, the integration constant may be fixed in terms of A and u. 

Clearly, the ^-coordinate corresponds to a non-compact extra dimension as the con- 
formal factor increases, or decreases, exponentially. The similarity with the extra fifth di- 
mension of the Randall-Sundrum model is striking: the size of this dimension can become 
finite only by introducing two branes (3-branes in the case of the RSI model, 4-branes in 
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this case) at two different points along the transverse ^-dimension. Then, the inter-brane 
distance defines the size of the extra dimension. The monotonic behaviour of the scale 
factor in the <fi direction, calls for the introduction of a pair of branes with positive and 
negative tension, as discussed earlier. We will show that this is indeed the case at the end 
of this section. Single-brane configurations could be also accommodated in our analysis by 
sending the second brane to an infinite distance from the first one. 

At this point, we can easily derive the relation between the internal 2D curvature and 
the bulk cosmological constant. The 2-dimensional curvature scalar is 

R(2) = Rl + R 6 6 = y ' ( 3 ' 12 ' ) 

From the solution ( |3.9|) , we further see that -j- = A 2 , so that 

R (2) = ^A B . (3.13) 

The behaviour of the warp factor along the ^-dimension, and subsequently the topology 
of this dimension, is strictly defined by the sign of the A 2 parameter which determines the 
sign of the two-space curvature through ( |3.12|) , or equivalently through eq. (|3.10|) , by the 



sign of the bulk cosmological constant. We now distinguish the two cases: 

(A): A 2 > 0. This case corresponds, through eq. (|3.10|) , to a negative bulk cosmological 



constant and to a negatively curved two-dimensional extra space-time 

A B <0, R {2) <0. (3.14) 



From the expression ( |3.9| ), we see that the function f(9) is characterized by the existence 
of a minimum at 

9 min = 9 ±j In 2a;, (3.15) 

under the assumption that f(9) is symmetric under the transformation 9 <-> —9. Identifying 
the two minima (for simplicity, we may set 9q — 0), we can compactify this extra dimension. 
Then, the quantity 2X9 min is the physical size of the extra dimension, 2(b 9 min ), in units 
of M 6 \ 

If we allow the two symmetric branches of the function f{9) to meet at 9 = 0, a cusp 
is inevitably created. One might be tempted to introduce an infinitely-thin 3-brane at 
this point creating a set-up that resembles the 5D single-brane configurations presented 
in Ref. M. However, the absence of any discontinuity in the first and second derivative 



8 



of the warp factor with respect to the ^-coordinate does not allow for the introduction of 
any infinitely-thin 3-brane in the model. The only allowed configuration is the one where 
each one of the 4-branes defining the size of the yj-dimension is a thick 3-brane extending 
along the ^-dimension similar to the one describing in Refs. [0, . In that case, the warp 
factor and its derivatives with respect to the ^-coordinate are everywhere well defined as 
demanded. Introducing a thick 3-brane along the ^-dimension will not spoil the solution 



( 3.11 ) found above: it would merely render it as the solution outside the "wall" of the 
3-brane where the minimum must take place. 

(B): A 2 < 0. This case arises under the assumption that the bulk cosmological constant 
takes a positive value. The internal curvature of the (<p, 6 l )-submanifold, in this case, is also 
positive according to the definition ( |3.12| ). Since A 2 < 0, we may write A = iX, where A 
is a real number. It turns out that for special values of the parameter u, the ^-dependent 
part of the solution for the warp factor ( |3.11D becomes periodic, and thus spontaneously 
compactified without the need for the introduction of any thin or thick 3-branes. However, 
in each case the resulting form of the line-element of the 6D spacetime contradicts basic 
assumptions of this analysis. 

Thus, if we choose u 2 = 1/4, the solution for f(9) is given in terms of a cos-type 
function. However, the presence of the coefficient 1/(4A) in eq. (|3.9|) renders this metric 
function imaginary leading to: 

ds 2 = a 2 {6, (p) {-dt 2 + dx 2 + dy 2 + dz 2 ) + b 2 [d6 2 - f{9) 2 d^ 2 } , (3.16) 

where we have set f(9) = if (9). As a result, the character of the yj-dimension changes from 
space-like to time-like which is in disagreement with our argument that the yj-dimension 
plays the role of the extra transverse dimension of the RS model. If, alternatively, we 
choose oj 2 = —1/4, the function f{9) comes out to be proportional to a sin-type function, 
however, in this case, u itself comes out to be imaginary. If we perform the following 
coordinate transformation 

tp^wi, (3 ' 17) 

the 6D line-element takes the form 

ds 2 = sin 2 [A (0 - 9 )} 1-dP + e 2Hi (dx 2 + dy 2 + dz 2 + ^L) j + 6 2 ^ 2 (3.18) 

where H = A/ (26 ) and where, for simplicity, we have set ao = 2A in eq. (|3.11|) . The above 
line-element describes a 6-dimensional spacetime whose 4 spatial dimensions are inflating. 
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This is clearly in contradiction with our main goal to derive solutions with static extra 
dimensions. 

We finally turn to the jump conditions that will help us determine the uj parameter of 
the solution as well as the number of fine-tunings that the model demands. Substituting 
the solution for ( |3.7|) into the jump conditions ([2.121) , we obtain the result 

±4^ = -Ai = A 2 - (3.19) 

The same fine-tuning between the two brane tensions, that was necessary for the RS solution 
to be consistent with the boundary conditions, appears also in our model. Note, however, 
that the second RS fine-tuning, between each brane tension and the bulk cosmological 
constant is absent. The value of A# defines the parameter A while uj is defined in terms of 
the brane tensions. No relation between A and uj exists in our model, and thus, Kb and Aj 
remain uncorrelated. Nevertheless, since only one of the two brane tensions has been fixed 
so far, another fundamental parameter of our model, if not the second brane tension or the 
bulk cosmological constant, should be fixed instead, in this case. Through eq. ( |3.15| ), the 
physical size of the ^-dimension is given in terms of the A and to parameters. By using eqs. 
( |3.15|) , and (|3.19|) , we obtain the result 

Qg^mi in 

^-2l- = l A ^ (3-20) 

where |Aj| stands for the absolute value of the brane tensions. The above relation fixes the 
physical size of the longitudinal extra dimension, which constitutes a fundamental parame- 
ter of the model, in terms of the brane tensions. The logarithmic dependence on the value 
of |Aj| ensures the smallness of the size of the longitudinal dimension even for large values 
of the brane tension. It is therefore the existence of this extra dimension that introduces 
an additional parameter in the model whose fixing replaces the fine-tuning between bulk 
and brane parameters. Note, however, that the locations of the two 4-branes along the 
transverse (^-dimension and, thus, the size of this dimension remain a free parameter. We 
address this point in the next section. 

A final comment on the number of fine-tunings is in order at this point: had the solutions 
for A < led to a consistent spontaneous compactification of the ^-dimension, an additional 
problem would appear: the special values of the uj parameter, for which these solutions 
arise, would result in the fixing of the values of both brane tensions through (|3.19|) . The 
fine-tuning between brane tensions and bulk cosmological constant would still be absent 
however the number of necessary fine-tunings in the model would be again two. We might 



10 



therefore conclude that periodic solutions, which respect the assumptions of our model and, 
at the same time, demand less fine-tuning than the usual 5D brane models, cannot arise in 
the framework of our analysis. 



4 Non-static 4D solutions 

Next, we proceed to construct 6-dimensional solutions with a 4-dimensional time-dependent 
submanifold but with a constant radion field once again. The source of this time-dependence 
will be a non- vanishing 4D effective cosmological constant. Solutions similar to these but 



in the presence of only one extra dimension have been derived before Here, we will 
investigate the possibility whether such solutions arise in the case of one additional, extra, 
longitudinal dimension. As in the previous section, we are going to assume that, initially, 
both extra dimensions are non-compact with the size of the yj-dimension becoming finite 
due to the presence of the two 4-branes while the ^-dimension will be appropriately warped 
and thus spontaneously compactified. 



Going back to the full time-dependent Einstein's equations (|2.3|) - (|2.9|) , we try once 
again the factorized ansatz: a(t,6,<p) = T(t) 0(#) The off-diagonal component (|2.9|) 

gives, as in the static case, a proportionality relation, 0(6) = Agf(8), between the two 6- 



dependent functions of the metric tensor. The remaining off-diagonal equations , eqs. (|2.7| )- 



(|2.8Q ) in conjunction with the above factorized ansatz, both lead to the result b = const., 
that guarantees the staticity of the extra 2-dimensional spacetime. 

Subtracting eq. fl2.6|) from eq. ( |2.5| ), we recover one of the two equations that determine 
the solutions for the metric functions $((p) and f(9), namely 

^ = ff"-f ,2 =u\ (4.21) 

Before trying to derive the second equation, we need to determine the solution for the 
time-dependent function T(t). Taking the sum of eqs. (|2.3| ) and (|2.4j) , we find the result 

T (~T~ T ' ~~ ^ (4.22) 

Obviously, one solution of the above equation is the trivial one, T = 0, which leads to 
the static case of the previous section with the Minkowski-like 4D submanifold. Clearly, 
the model accepts another solution that may be determined by demanding the expression 



11 



inside brackets to be zero. Then, we obtain the alternative solution 

where cq and to are integration constants. If we pass from the conformal time t to the 
physical time t, through the transformation dt = T(t) dt, we might easily see that the 
solution takes the form 

Tit) = e H{i ~ io) . (4.24) 

In the above expression, H = cq is again an integration constant, which may be either real 
or imaginary describing a de Sitter 

ds\ = -dt 2 + e 2Hi (dx 2 + dy 2 + dz 2 ) (4.25) 

or Anti de Sitter 

ds\ = dx 2 + e 2Hx (-dt 2 + dy 2 + dz 2 ) (4.26) 

4D submanifold, respectively. Since the Einstein's equations are all expressed in terms of 
the conformal time, we will continue using t, instead of t, and the form of the solution ( |4.23|) , 



instead of Q4.24 ), for convenience. In terms of the conformal time, a real or imaginary cq 



will distinguish between a de Sitter or Anti de Sitter 4D submanifold. 

Having determined the solution for T(t), we now turn to eq. fl2.6Q , which can be brought 



to the form 



Comparing the above with eq. (|3.6|) , we may easily see that we can recover the equation 



^TS| ) for f(9), if satisfies the following equation 



The solution for the (^-dependent part of the warp factor, now, takes the form 

$(^) = -^ sinh M^-^°l]' (4.29) 

for de Sitter space-time, and 

= cosh[w ^ _ ^ o) ] (4.30) 

for Anti de Sitter space-time. 
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Since the basic equation for the metric function f{9) has remained unchanged, the 
general solution given by eqs. (|3.9|) and (|3.10|) still remains the same. The subsequent 
discussion on the possible ways of compactifying the ^-dimension, in the case of negative or 
positive bulk cosmo logical constant, and the main conclusions drawn at the end of section 3 
still hold. On the other hand, the time-dependence of the warp factor has radically changed 
the solution for its (^-dependent part. Instead of the monotonic exponential dependence 
that prevailed in the static case, the time-dependent case may accommodate both sinh and 
cosh-type solutions. Solutions similar to the above usually arise in the case where classical 
19| or quantum |20| effects from bulk scalar fields are taken into account or when the effect 



of a bulk stabilizing potential is included in the model []B5[ ^ . Note, however, that no such 
effects have been assumed to be present in our analysis. The conclusions drawn in the 
previous section from the jump conditions are also likely to change. Substituting the above 
solution for $(<£>) into the jump conditions fl2.12| ), we obtain the constraints 

k 2 

cu coth[w fa - <p )] = (-I) 1 b A, (4.31) 
for the sinh-type solutions, and 

K 2 

u tanh[u; - <p )] = (-1)* -f b A, (4.32) 

for the cosh- type solutions. In the above two equations, 2=1,2 denotes again the two 
4-branes. The sinh-type solution with its monotonic behaviour can clearly accommodate 
only pairs of positive-negative tension branes in analogy with the static case. Moreover, 
the specific solution is plagued by the existence of a singularity at the point tp — (fo- In 
order to exclude the singularity from the 6D spacetime, both branes must be located on the 
same size of the singularity, i.e. ipi,ipz < y^o- On the other hand, the cosh-type solution is 
everywhere well defined. The same point, ip = <y9 , corresponds, in this case, to a minimum 
which allows for pairs of positive tension branes to fit in the model. In both cases, either 
for sink or cosh-type solutions, the above jump conditions will fix the location of the two 
branes relative to the singularity or the minimum, respectively, in terms of the two brane 
tensions. Note, however, that the brane tensions, considered as input parameters of the 
model, remain totally uncorrelated. Moreover, the lack of fine-tuning between any of the 
brane tensions and the bulk cosmological constant still holds rendering the model free of 
any fine-tuning. 
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5 Radion dynamics 



The main goal of the previous sections was to find 6D solutions, static and non-static, with 
a constant radion field which demand less or no-fine-tuning of their parameters compared to 
other models in the literature. Nevertheless, an important aspect of these solutions needs 
to be studied next: do these solutions actually extremize the radion effective potential? 
And, if yes, is this extremum a minimum or a maximum of the effective potential? 



In Ref. [IIJ , a 5-dimensional "extremization" constraint that could serve as a consistency 
check for any solutions with a constant radion field, was derived. It might be worth deriving 
the corresponding constraint in 6 dimensions and comment on the possible differences that 
arise as one changes the number of extra dimensions. In order to do that, we need to go 
back to the time- dependent Einstein's equations, and start by constructing the 4D trace of 
the energy- momentum tensor (by taking the sum of eq. ( |2.3|) with three times eq. (|2.4| )), 
which comes out to have the form 



4th 



46 2 



a 2 b 2 



6b \2ab 6a 4/'' 
a 2 b a 3 b a 3 

12 d 2 a 



+ 



12 



+ 



+ 



b 2 f b 2 f 2 a 2 
12 f da 12 



da\ 2 
d^J 



da\ 2 12 d 2 a 



+ 



(5.33) 



b 2 f 2 a dcp 2 b 2 fa d6 b 2 a 2 \d9 J b 2 a d6 2 ' 
Next, we construct the following linear combination of all the components of the 6D energy- 
momentum tensor 



h i 



(4-n)r; + 3(n-2)T 5 5 + 3(n-2)T 6 f 



. 46 2 126 24d& 
a z ¥ a l b a 4 b 



12(n-l)- + (4 

a 6 



n) 



.4T 
b 2 f 



+ 



24 



[n 



1) (^)\ J 

f 2 a 2 \difj f 2 a dip 2 



d 2 a f da 

+ h 

fa 09 



n 



da\ 2 



1 d 2 a 



(5.34) 



a" \ou/ a d6 2 . 

where n is an integer. In the case where the warp factor can be written as the product of 
two functions, one depending on the 4D coordinates and one on the extra coordinates, the 
above expression can be greatly simplified. If we write the warp factor in the factorized 
form 

a(t, 9, (p) = T(t) W(9, (p) , (5.35) 
the expression in the last line of eq. ( |5.34j ) may be written as 



24 



(n — 1) (da 



Pa 2 



dip) f 2 a dip 2 



d 2 a f da 

+ h 

fad9 



1) /da\ 2 1 d 2 a 



(n - 1) /da\ 



+ 



deJ a de 2 
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24 



1 



nf x m ) 



24 



(D m D m W r 



(5.36) 



where x m denotes the extra coordinates (9, if) and g^ n the metric tensor of the 2-dimensional 
extra spacetime. The above 2D double covariant derivative of the function W n divided by 
W n is the analogue of the double derivative of the same function appearing in the 5D 
constraint of Ref. flHj. In the more general case where more than one, possibly non-flat, 
extra dimensions are present in the theory, the double derivative needs to be replaced with 
the covariant derivative in order for the internal geometry of the extra spacetime to be 
taken into account. 



In addition, we may rewrite some of the other terms appearing in eq. ( 5.34 ) in the 
following way 

a 2 4/" 

where Ru\ and Rm stand for the 4D and 2D, respectively, scalar curvature. Then, the 
constraint (|5.34j) takes the simplified form 



24 



(D m D m W r 



(4 



3(n - 2) Tr + 3(n - 2) T fi fc 



+ (n-l)— R (A) + 2(4 -n)R (2) + 



(4 



n 



4b 2 12b 
a 2 b 2 a 2 b 



24ab 



(5.38) 



Compared to the 5-dimensional one |10[, the 6-dimensional version of the above constraint 
has a similar but more generalized structure. It involves all the extra components of the 
energy- momentum tensor, namely T 5 5 and T| , as anticipated, and both on an equal footing. 
Moreover, the 2-dimensional scalar curvature of the extra spacetime explicitly makes its 
appearance together with the 4-dimensional one. Finally, the coefficients appearing in 
front of the 4D trace and extra components of the energy-momentum tensor seem to be 
"dimension" -dependent. By mere comparison of the 5D and 6D versions of the constraint, 
we may easily conclude that the coefficient in front of the extra components behaves as 
(1 + d)(n — 2), where d the number of extra dimensions. However, the dependence of the 
coefficient of the 4D trace is more subtle and a higher dimensional calculation could only 
reveal its exact form in terms of d. 

We would also like to stress here an additional point on the form of the constraint ( 5.38|) . 
It holds for every solution of the 6D Einstein's equations, even for the ones with a non-static 
extra spacetime. For the particular case of n = 1 the above constraint can be interpreted 
as the equation of motion of the time-dependent scalar field b(t). This can become clear if 
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we rewrite Eq. ( |5.38|) as 




b 2 = 4 (T; - if - T 6 6 ) + 2R (2) 



8 

w 



(D m D m W) . 



(5.39) 



For solutions with a non-static extra spacetime, the rhs of the above equation vanishes 
when an extremum, either minimum or maximum, is reached. For solutions with a constant 
radion field, which by definition correspond to an extremum of the radion effective potential, 
the same combination should also vanish. Therefore, the solutions found in the two previous 
sections, either static or non-static, but with constant b, should satisfy the above constraint. 
By using the following expressions for the components of the energy-momentum tensor 



T; = -4A B -— T 5 5 = -A B - — Lsfr-ipt), Ti = -A B , (5.40) 



setting W = 0(6) *&((p) and employing the form of the solutions found in sections 3 and 4, 



By using the jump conditions ( |2.12| ) and the definition of the A parameter from eq. ( |3.10| ), 
we may easily conclude that both the static and non-static solutions derived in the previous 
sections satisfy the "extremization" constraint, as anticipated. 

In order to answer the question whether the above extremum is a minimum or a maxi- 
mum, a perturbation analysis around the above solutions needs to be performed, in which 
the time-dependent, small perturbation will be associated with the radion field. Therefore, 
we consider the following ansatz for the line-element of the 6-dimensional spacetime 



Note that we have switched to the system of non-conformal coordinates and we have set 
b = 1, for simplicity. We have also chosen to perturb the de Sitter solutions found in the 
previous section, however, our analysis can be easily extended to the cases of Minkowski or 
Anti de Sitter branes by setting if 2 — > or H 2 — > —H 2 (with the t-dependent perturbations 
replaced by x-dependent), respectively. 

The above ansatz when substituted into the 6D Einstein's equations will lead to a 
system of differential equations and constraints that govern the behavior of the new line- 
element. We will work in the linear order approximation and, thus, we are keeping only 



we obtain 





ds 2 = [a 2 f(6) 2 $(<p) 2 + e A(6, if) b(t)} (-dt 2 + e 2Ht dx 2 ) 




(5.42) 
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terms proportional to the small parameter e< 1. In this approximation, the generalized 
(00), (ii), (05) and (06)-components of Einstein's equations take the form 

B + ^p) *Hb - (4alf^A B - 3*) Bb-^f-^ 

3bf'dA 3bd 2 A a 2 b$ 2 fd 2 B dB , f 2&B\ _ 2 

B + -^-) (-2Hb -b) + (4alf^A B - 3H 2 ) Bb + ^f- + 36 ^ 1 



alf 2 <& 2 ) v ; v 6 0J " ' f 2f 2 dip 2 

3bf OA 3bd 2 A a 2 b<& 2 (d 2 B ££l dB , 2 d 2 B\ 2 . . AA , 

"Ww + ^ + ^W + ff w + f V) = ~ KlAh Ab ' (5 - 44) 

3A \f 3 dA 1 dB _ 

3A \ $' 3 &4 1 &B . , 

a 2p$2j $ 2a( ) f 2 § 2 dip 2 d<p 

respectively. 

By taking the sum of the (00) and (^)-components, we end up with the constraint 

( fl + s^p)( jri - s ) = °- < 5 - 47 > 

which demands that one of the two expressions inside brackets vanishes. If we assume that 
the expression inside the first bracket is zero, then, we can determine the exact form of the 
unknown functions A(8, ip) and B(9, ip) by plugging this constraint into the off-diagonal 
components ( |5.45|) - (p.46|) . Then, we find that 

A V ,<„) = 1, B(B, V ) = - ^ (5.48) 



modulo a multiplicative constant. Both of the equations Q5.43D - (|5.44|) , in that case, reduce 



to the background equation Q4.28|) (with Cq being replaced by H 2 due to the use of non- 
conformal coordinates) which is obviously satisfied by the background solution. However, 
the solution ( |5.48|) fails to satisfy any of the remaining components of Einstein's equations. 
For example, the off-diagonal (56)-component, which has the form 

2f (OA &\ (f dB & dB^ & dA d 2 A _ 

T W " o/ \1 + ¥ -bJ) + ¥ m ~ defy ~ ' (5 ' 49) 
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leads to the constraint /'$' = which is obviously in contradiction with the form of the 
background solution. Finally the (55) and (66)-components, which after using eq. (|5.48j) 
take the simplified form 

(D^b) I2b(& 2 H 2 \ 12bf' 2 46$" 2 

' 1 -K e bA B , (5.50) 



a §/2$2 p y $ 2 fl 2 $2 ; p p§ 

{Djjyty b (W 2 \2H 2 \ Ubf 2 46/" _ 2 
a 2 f 2 <S> 2 P\W "a 2 ® 1 ) — - K « fiAfl ' V- bL > 

respectively, also turn out to lead to an inconsistent result. As a result, the option that the 
first expression inside brackets in eq. (|5.47|) is zero needs to be excluded. We are therefore 
left with the alternative option (Hb — b) = which leads to an exponential solution of the 
time-dependent perturbation. By employing the equation of motion of the "radion" field 

g^D^DJ) = m 2 b (5.52) 

and the exponential form of its solution, we get the result m 2 = —AH 2 . Obviously, the 
mass squared of the radion field turns out to be negative for de Sitter 4D subspace, zero 
for Minkowski and positive for Anti de Sitter 4D subspace. 

The above result is in perfect agreement with similar works conducted in 5 dimensions: 
in Refs. [[|, [Tl| (see also |12|), it has been shown that the system of two Minkowski 



branes with a zero total brane tension cannot be stabilized since it leads to a radion field 
with a zero mass. In that case, the extremization of the radion potential follows from the 
fine-tuning of the parameters of the model which however cannot create a unique minimum 
for the radion field in the absence of a physical stabilization mechanism. As a result, the 
static solution derived in Section 2 describing a pair of flat branes cannot be stabilized and 
the radion field remains always massless. The same conclusion was drawn in Refs. |I3L |I4 



where the stability of curved branes was also studied. According to their results, a pair of de 
Sitter branes leads to an effective theory for the radion field with a negative mass squared 
while a pair of Anti de Sitter branes turns out to be stable since it leads to a positive radion 
mass squared. Therefore, our solution (|4.29|) corresponding to two branes with a positive 
effective cosmological constant is unstable under small time-dependent perturbations while 
the alternative solution described by ( |4 . 3 Q|) and corresponding to two branes with negative 



effective cosmological constant is stable. We bear in mind that this conclusion ignores 
any additional contribution to the 4D radion potential as might be expected from a more 
complete model which includes supersymmetry and supersymmetry breaking. Let us finally 
note that the above results are also in agreement with those of Ref. [0] where the extremum 
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of the radion effective potential was a global minimum only in the case of a 4D negative 
cosmological constant. 

6 Conclusions 

In this paper, we have presented a 6D brane world model in the framework of which we 
addressed a number of important issues that arise in the context of higher-dimensional brane 
models. The relaxation of the fine-tuning of the fundamental parameters of the model was 
one of them: solutions with less or no fine-tuning at all, compared to their 5D analogues, 
were constructed where the severe correlation between brane and bulk parameters was 
replaced by the fixing of the sizes of the two extra dimensions. In this way, another problem, 
that of the determination of the inter-brane distance was automatically resolved. The 
solutions derived allow for the introduction of pairs with positive-negative brane tensions, 
in analogy with the RS model, however more physically interesting configurations with 
pairs of only positive tension branes were also found. Finally, the issue of the behaviour of 
the aforementioned solutions, under time- dependent perturbations around configurations 
with a constant radion field, was examined and conclusions regarding their stability were 
derived. 

In more detail, considering a 2D internal space of constant curvature, an exact static 
solution was first presented. The warp factor exhibits neither spherical nor cylindrical 
symmetry but depends on both extra coordinates. Two 4-branes are introduced at two 
different points along the so-called "transverse" extra dimension, along which the warp 
factor is a pure exponential resembling the profile of the warp factor in the case of the 
5D RS model. The remaining extra dimension, the "longitudinal" one, along which the 
4-branes extend, was shown to be "spontaneously" compactified in the case of a negative 
bulk cosmological constant, or equivalently of a negatively curved internal space. Although 
the consistency of the bulk solution with the brane boundary conditions demand the two 
branes to have exactly equal and opposite tensions, no correlation exists between the brane 
tensions and the bulk cosmological constant. Instead, it is the size of the longitudinal extra 
dimension that is fixed through the jump conditions in terms of the value of the brane 
tensions. 

The above solution, although it is characterized by reduced fine-tuning of its funda- 
mental parameters, has the inter-brane distance along the transverse extra dimension as a 
free parameter. In an attempt to resolve this problem, too, we then derived a non-static 
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solution whose 4D subspace corresponds to a de Sitter or Anti de Sitter spacetime. These 
solutions have a number of positive features: first, being non-static due to the presence of 
an effective 4D cosmo logical constant, the only fine-tuning of the model, the one between 
the two brane tensions, is also relaxed rendering the model free of any fine-tuning; sec- 
ond, although they have the same profile along the compactified, longitudinal dimension, 
the exponential behaviour along the transverse one changes to a cosh or sinh-tjpe one, a 
fact which has two important consequences: the fixing of the locations of the two branes 
and, thus, of the inter-brane distance, and the accommodation of pairs of positive tensions 
branes instead of only pairs of positive-negative tensions. 

The fixing of the physical inter-brane distance and thus of the size of the extra dimen- 
sion relies on the assumption that the scale factor along the internal spacetime remains 
constant. Both of the above solutions, static or non-static, were derived under this assump- 
tion, or, alternatively, that the corresponding "radion" field was already at the extremum 
of its effective potential. An "extremization" constraint, valid for solutions with static or 
non-static radion field whose effective potential possess an extremum, was formulated and 
used as a consistency check for our solutions. The final issue of whether this extremum 
was a unique minimum or merely a local maximum of the radion effective potential needed 
to be addressed. The background solutions with a constant scale factor along the extra 
dimensions were perturbed and a linear stability analysis was performed. Our analysis 
revealed the stability of the non-static solutions describing an Anti de Sitter 4D spacetime 
and accommodating pairs of only positive tension branes. The remaining solutions describ- 
ing 4D de Sitter and Minkowski spacetimes were found to correspond to local maxima and 
saddle points, respectively, in close analogy to similar analyses performed in 5 dimensions. 
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